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Formulas for calculating the extremal ranks and inertias of a 
matrix- valued function subject to matrix equation restrictions 
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Abstract. Matrix rank and inertia optimization problems are a class of discontinuous optimization problems in which 
the decision variables are matrices running over certain matrix sets, while the ranks and inertias of the variable matrices 
are taken as integer-valued objective functions. In this paper, we establish a group of explicit formulas for calculating the 
maximal and minimal values of the rank and inertia objective functions of the Hermitian matrix expression Ai — BiXB* 
subject to the common Hermitian solution of a pair of consistent matrix equations B2XB2 = A2 and B3XB3 — A3, 
and Hermitian solution of the consistent matrix equation B4X = A4, respectively. Many consequences are obtained, in 
particular, necessary and sufficient conditions are established for the triple matrix equations BiXB* — A\, B2XB2 = A2 
and B3XB3 = A3 to have a common Hermitian solution, as necessary and sufficient conditions for the two matrix equa- 
tions BiXB* — Ai and B4X — A4 to have a common Hermitian solution. 
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1 Introduction 

Throughout this paper, 



q 

<-j ' £mxn s t anc j s f or th e set of all to x n complex matrices; 

^ 1 stands for the set of all m x m complex Hermitian matrices; 

A T , A*, r(A), M(A) and ^V(A) stand for the transpose, conjugate transpose, rank, range (column space) 
and null space of a matrix A G C mx ™ 5 respectively; 

I m denotes the identity matrix of order to; 

[A, B] denotes a row block matrix consisting of A and B; 

A > (A ^ 0) means that A is Hermitian positive definite (Hermitian positive semi-definite); 



two A, B G are said to satisfy the inequality A > B (A ^ B) in the Lowner partial ordering if A — B 
is positive definite (positive semi-definite) ; 

the Moore-Penrose inverse of A € C mxi \ denoted by A\ is denned to be the unique solution X satisfying 
the four matrix equations AXA = A, XAX = X, (AX)* = AX and (XA)* = XA, which satisfies 
AAt = A^Aii A = A*; 

a matrix X is called a Hermitian g- inverse of A G Cjj, denoted by A~ , if it satisfies both AXA — A and 

X = X*; 

Ea and Fa stand for Ea = I m ~ AA^ and Fa = I n — A^ A. The ranks of Ea and Fa are given by 
t(Ea) = to — r(A) and t(Fa) = n — r(A); 

i+(A) and i- (A), usually called the partial inertia of A G Cjj, are defined to be the numbers of the positive 
and negative eigenvalues of A counted with multiplicities, respectively, which satisfy r(A) = i + (A)+i_ (A). 

The matrix approximation problem is to approximate optimally, with respect to some criteria, a matrix by 
one of the same dimension from a given feasible matrix set. Assume that A is a matrix to be approximated. 
Then a conventional statement of general matrix optimization problems of A from this point of view can be 
written as 

minimize p( A — X ) subject to X G S, 0-A) 

where p(-) is certain objective function, which is usually taken as the determinant, trace, norms, rank, inertia 
of matrix, and S is a given feasible matrix set. A best-known case of is to minimize the norm || A — X |||, 
subject to X G S. 
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In this paper, we take the matrix set S as 

S = {4>{X) = A 1 -B 1 XBl | [B 2 XB* 2 , B 3 XB%] = [A 2 , A 3 ]}, (1.2) 
S = { 4>(X) = Ax- B t XB* | B 4 X = A4, }, (1.3) 

where A t G and B 4 G C" uxn , Aa, 5 4 G C mx ™ are given, i = 1, 2, 3, and X € C£ is a variable matrix, and 
study the following constrained optimization problems: 



Problem 1.1 For the constrained linear matrix- valued function in (|1.2[) . establish explicit formulas for calcu- 
lating 

maxr(A 1 -B 1 XBl) s.t. [B 2 XB*, B 3 XB*] = [A 2 , A 3 ], (1.4) 

min r(A 1 -B 1 XBf) s.t. [B 2 XB* 2 , B 3 XB* } = [ A 2 , A 3 ], (1.5) 

^ H 

mm i ± {A l -B 1 XB* 1 ) s.t. [5 2 XB*, B 3 XS| ] = [ A 2 , A 3 ], (1.6) 

H 

mmi^Ax-BtXBl) s.t. [MB^, B 3 XB* 3 } = [A 2 , A a ]. (1.7) 



Problem 1.2 For the constrained linear matrix- valued function in (|1.2[) . 

(i) establish necessary and sufficient conditions for the following three matrix equations 

[BiXBl B 2 XB* 2 , B 3 XB* 3 ] = [A u A 2 , A 3 ] (1.8) 
to have a common Hermitian solution; 

(ii) establish necessary and sufficient conditions for A\ — B\XB\ > (^, <, ^)0 to hold for an X 6 
satisfying [ B 2 XB* 2 , B 3 XB* 3 ] = [ A 2 , A 3 ]; 

(iii) establish necessary and sufficient conditions for A\ — B\XB\ > (^, <, ^)0 to hold for all X G 
satisfying [B 2 XB 2 , B 3 XB 3 ] = [A 2 , A 3 ], namely, A\ — B\XB\ is a positive map under the restriction 
[B 2 XB* 2 , B 3 XB%} = [ A 2 ,A 3 }. 



Problem 1.3 For the constrained linear matrix- valued function in (jT73j) , establish explicit formulas for calcu- 
lating 

maxrMi -B^B*,) s.t. B A X = Aa, (1.9) 

min r(Ai - B X XBX ) s.t. B 4 X = Aa, (1.10) 

xecv. 



s.t. 


B A X 


= Aa 


s.t. 


BaX 


= Aa 


s.t. 


BiX 


= Aa., 


s.t. 


BiX 


= A 4 . 



maxi±(Ai ~B x XB;) s.t. BaX = Aa. (1.11) 
mm i±( A Y - B X XB* ) s.t. BaX = Aa. (1.12) 

Problem 1.4 For the linear matrix map in (| 1 . 3[) . 

(i) establish necessary and sufficient conditions for the following two matrix equations [BiXB^, B4X] — 
[A±, A4] to have a common Hermitian solution and nonnegative definite solution; 

(ii) establish necessary and sufficient conditions for A\ — B\XB\ > (^, <, ^)0 to hold for an X G 
satisfying B4X — A4; 

(iii) establish necessary and sufficient conditions for A\ — B\XB\ > (^, <, ^)0 to hold for all X G 
satisfying B4X = A4, namely, Ai — B\XB* is a positive map under the restriction B4X = A4. 

The extremal ranks and inertias of a matrix expression can directly be used to describe some behaviors of 
the matrix expression, for example, 

(I) the maximal and minimal dimensions of the row and column spaces of the matrix expression; 

(II) nonsingularity of the matrix expression when it is square; 

(III) solvability of the corresponding matrix equation; 

(IV) rank, inertia and range invariance of the matrix expression; 

(V) semi-defmitcness of a matrix expression; etc. 

On the other hand, matrix rank and inertia optimization problems are NP-hard in general due to the discon- 
tinuity and combinational nature of rank and inertia of a matrix and the complexity of algebraic structure of 
S. 
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Mappings between matrix spaces with symmetric patterns can be constructed arbitrarily, but the linear map 
4>{X) in (|1.2[) and (|1.3|) is the simplest cases among all LMFs with symmetric patterns. This matrix-valued 
function is the starting point in dealing with various complicated matrix-valued functions with symmetric 
patterns. In recent years, the matrix-valued function </>(X) = A — BXB* was reconsidered and many new 
results on its algebraic properties were obtained, for instance, 

(i) Expansion formulas for calculating the (global extremal) rank and inertia of <fi{X) when X running over 
Cg, see [TH [301 [38]. 

(ii) Nonsingularity, positive definiteness, rank and inertia invariance, etc., of the 4>(X), see [30l 138) . 

(iii) Canonical forms of the <f>(X) under generalized singular value decompositions and their algebraic proper- 
ties, see [19]. 

(iv) Solutions and least-squares solutions of the matrix equation <j>(X) — and their algebraic properties, see 

US EH 1331031371. 

(v) Solutions of the matrix inequalities 4>(X) > (^, <, and their properties, see [30] , 

(vi) Minimization of tr[0(X)0*(X) ] s.t. r[4>(X)] = min, see [37] , 

(vii) Formulas for calculating the extremal rank and inertia of <j>{X) under the restrictions r{X) ^ k and/or 
±X > 0, see [35]. 

(viii) Formulas for calculating the extremal rank and inertia of cf>(X) subject to a consistent matrix equation 
CXC* = D, see [E]. 

This basic work was also extended to some general LMFs, such as, A — BX - (BX)*, A - BXB* - CYC* and 
A — BXC — (BXC)* , where X and Y are (Hermitian) variable matrices of appropriate sizes; see El [T51 [HJ IT71 

HH1 133 1331 ESI- 

We shall use some pure algebraic operations on matrices to derive two groups of analytical formulas for 
calculating the global extremal values of the objective functions in (|1.4j) - (|1.7j) and (|1.9|) - (|1.12|) . and then to 
present a variety of valuable consequences of these formulas. 

Since variable entries in a matrix-valued function are often regarded as continuous variables in some con- 
strained sets, while the objective functions — the rank and inertia of the matrix- valued function take values only 
from a finite set of nonnegative integers, Hence, (ll.4|) - (jl.7|) and (|1.9[) - (|1.12[) can be regarded as continuous- 
integer optimization problems subject to equality constraints. This kind of non-smooth optimization problems 
cannot be solved by using various optimization methods for solving continuous or discrete cases. There is no rig- 
orous mathematical theory for solving a general rank and inertia optimization problem due to the discontinuity 
and nonconvexity of rank and inertia of matrix. In fact, it has been realized that rank and inertia optimization 
problems have deep connections with computational complexity, and are regarded as NP-hard in general; see, 
e.g., [II El 111 [5l [3 111 [9l [12l [22 [25l [27] . Fortunately, some special rank and inertia optimization problems now can 
be solved by pure algebraical methods. In particular, analytical solutions to the rank and inertia optimization 
problems of the 4>(X) in (fL~2|) and (fL~3]> . as well as (fTTl jl — (fTTT j) and (fl~9]) - (|1.12|) can be derived algebraically by 
using generalized inverses of matrices. 

The following are some known results for ranks and inertias of matrices and their usefulness, which will be 
used in the latter part of this paper. 

Lemma 1.5 (|30j) Let H be a matrix set in . Then, 

(a) T-L has a matrix X > (X < 0) if and only if maxxefi i+(X) = m (maxxg^ i^(X) = m) . 

(b) All l£M satisfy X > (X < 0), namely, % is a subset of the cone of positive definite matrices (negative 
definite matrices), if and only if minx^-H i+{X) = m (minxey. i— {X) — m ) ■ 

(c) H has a matrix X ^ (X ^ 0) if and only if minxeH i—(X) = (minxe-H = 0) . 

(d) All X G H satisfy X (X ^ 0) namely, % is a subset of the cone of nonnegative definite matrices 
(semi- definite matrices), if and only if maxxeu i-(X) = (niaxxeH *+( X) = 0) . 

The question of whether a given function is (definite or semi-definite everywhere is ubiquitous in mathematics 
and applications. Lemma ll.5f aWd) show that if some explicit formulas for calculating the global maximal and 
minimal inertias of a given Hermitian matrix- valued function are established, we can use them, as demonstrated 
in Sections below, to derive necessary and sufficient conditions for the Hermitian matrix-valued function to be 
definite or semi-definite. 
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Lemma 1.6 Q21J) Let A E 

hold 



B E C mxp and C E C 9X ". Then, the following rank expansion formulas 



r[A, B] = r(A) + r(E A B) = r(B) + r(E B A), 

'A' 

C 

A B 
C 



= r(A) + r{CF A ) = r(C) + r(AF c ), 
= r{B)+r(C)+r{E B AF c ). 



Three useful rank expansion formulas derived from (|1.15|) are 

r(P) + r 



ABO 
COP 

A B ' 
C 
Q 

ABO' 
COP 
Q 



A B 
E P C 



A BFq 
C 



r(Q) + r 
r(P)+r(Q)+r 



A BF Q 
E P C 



We shall use them in Section 2 to simplify ranks of block matrices involving Ep and Fq. 
Lemma 1.7 ( [30j ) Let A E C£\ B E C mxn , D E Cg, and let 

U = 



' A 


B 


, v = 


' A 


B' 


B* 





B* 


D 



Then, the following expansion formulas hold 

i±(U)=r(B)+i±(E B AE B ), 

E A B 



(a) If A^fO, then 

(b) // A < 0, then 



i±(V) = i±(A) +i± b*Ea D-B*A^B 
i+{U) = r[A, B], i-(U)=r(B), r(U) — r[A, B] + r(B). 



i+(U) = r(B), i-(U) = r[A, B], r(U) — r[A, B] + r(B). 

(c) IfM(B) C M(A), then 

i±{V) = i±(A)+i±(D- B*A^B), r(V) = r{A)+r{D- B*A^B). 

(d) IfM{B)f\M{A) = {0} andM{B*)nM{D) = {0}, then 

»±(V) = i±(A) + i ± {D) + r(B), r(V) = r(A) + 2r(B) + r(D). 

Three general expansion formulas derived from (|1.19|) are 



A BF P 
F P B* 



= i± 



ABO 
B* P* 
P 



-r(P), 



A BF P 
F P B* 



ABO 
B* P* 
P 



We shall use them to simplify the inertias of block Hermitian matrices that involve Fp = I — P'P. 
Lemma 1.8 Let A 3 E C m J xn , B 3 E C pxq > and C 3 E C"^ x ^ be given, j = 1,2. Then, 
(a) |26j The pair of matrix equations 

AiXBx = Ci and A 2 XB 2 = C 2 
have a common solution for X E C™ xp if and only if 



3t{pi) C M(Aj), 3%(C*) C M(B*), r 








At 





-c 2 


A 2 




B 2 






A x 

An 



(1.13) 
(1.14) 

(1.15) 



(1.16) 
(1.17) 

(1.18) 



(1.19) 
(1.20) 

(1.21) 
(1.22) 
(1.23) 
(1.24) 



-2r(P). (1.25) 



(1.26) 



r[B u B 2 ), j = l,2. (1.27) 
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(b) |29) Under (|1.45|) . the general common solution to (jl.44l) can be written in the following parametric form 

X = X + F A V 1 + V 2 E B + F Al V 3 E B2 + F A2 V 4 E Bx , (1.28) 



where A 



At 
A 2 



B = [Bi, B2], and the four matrices Vi, . . . , V4 G C™xp are arbitrary. 



Lemma 1.9 Let A G C mxn and B G Cg- be given. Then, 

(a) [SUn] The matrix equation AX A* — B has a solution X G if and only if £%{B) C M{A), or equivalently, 
AA^B = B. 

(b) |30) Under AA^B = B, the general Hermitian solution of AX A* = B can be written in the following two 
forms 

X = A*B{A*)* + U-A*AUtfA, (1.29) 
X = A*B(A*)* + F A V + V*F A , (1.30) 

where U G and V G C nx ™ are arbitrary. 

More results on properties of solutions of AX A* = B can be found in [16l [20] . 

Lemma 1.10 ([10]) Let A, B G C mxn 6e given. Then, 

(a) TTie matrix equation AX = B has a Hermitian solution X G if and only if £tf,(B) C Sf\A) and 
AB* = BA*. Ln this case, the general Hermitian solution of AX = B can be written as 

X = A^B + (A^B)* - A^BA^A + F A UF A , (1.31) 

where U G C§ is arbitrary. 

(b) The matrix equation AX — B has a solution X G if and only if S%{B) C ,5? (A), A_B* and 
r(AB*) = r(£?). In this case, the general solution of AX — B can be written as 

X = B*(AB*) f B + F A UF A , (1.32) 

where U G Cg is arbitrary. 
Lemma 1.11 Le£ A G and B e C mxn be given. Then, 
(a) |30U38) The global maximal and minimal ranks and inertias of A — BXB* subject to X G Cj^ are given by 

■[A, B], (1.33) 



max r( A — 


BXB*) 


min r( A — 


BXB*) 


max i+( A — 


BXB*) 


min i± ( A — 


BXB*) 



A B 
B* 



A B 
B* 



A B 
B* 



(1.34) 
(1.35) 
(1.36) 



(b) [35] The global maximal and minimal ranks and inertias of A — BXB* subject to ^ X G C§ are given 
by 

(1.37) 
(1.38) 
(1.39) 
(1.40) 
(1.41) 
(1.42) 



max r( A 

o<xecg 


BXB* ) = 


--r[A, B], 


min r( A + 

0<xecg 


BXB* ) = 


i+{A)+r 


A, B]-U 


-(M), 


max i + ( A 


+ BXB* ) 


= i+(M), 


min i + ( A - 


{- BXB* ) = 


--i+(A), 






max i_( A 


+ BXB* ) 


= i-(A), 


min i-( A - 


\- BXB* ) -- 


= r[A,B}- 


-i+(M), 




max r( A - 


BXB* ) = 


--r[A, B], 


min r( A — 

0<xecg 


BXB* ) = 


+ r 


A, B] - L 


(M), 


max i+( A 


- BXB* ) 


= i+(A), 


min i + ( A - 


- BXB* ) = 


= r[A, B] - 


-»_(M), 




max i—( A 


- BXB* ) 


= i-(M), 


min i-(A 


- BXB* ) 


= i-(A). 
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Lemma 1.12 ([18j) Let AeC%, B £ C mxp and C G C qxm be given, and let 



M 1 = 



A B 
B* 



M, = 



'A 


C* ' 














n 
u 


5 










" A 


B 


c*~ 


, N 2 = 


'A 


B 


C*~ 


B* 








C 









N = [A, B, C* ], iVi 

Then, the global maximal and minimal ranks and partial inertias of A — BXC — (BXC)* are given by 
max r[A — BXC - (BXC)* 1 = min{r(N), r(i\M, r(N 2 )}, 

min r[A- BXC - (BXC)*} =2r(N) +max{s u s 2 , s 3 , s 4 }, 

max i±[ A - BXC - (BXC)* } = min{i±(Mi), i±(M 2 )}, 



min j±[ A - BIG - (MC)* 1 = r(iV) + max{ i±(Mi) - r(JVA i±(M 2 ) - r(jV 2 ) }, 



(1.43) 
(1.44) 

(1.45) 
(1.46) 
(1.47) 
(1.48) 



si = r(Mi) - 2r(JVi), s 2 = r(M 2 ) - 2r(iV 2 ), 
s 3 = i + (M x ) + t _(M 2 ) - r(JVi) - r(N 2 ), s 4 = ^(Mx) + i_(M 2 ) - r(JV x ) - r(JV 2 ). 
In particular, ifM(C*) C M(B), then 



max r[ A - B1C - (BXC)* 
min r[4- B1C - (BXC)* 



max i±[ A - B1C - (BXC)* 

xeCf x « 



min i±f A - B1C - (BXC)* 



min ^?'[ A, _B ], r 
i± 

r[A, B]+i± 



A C* 
C 

A C* 
C 



2r 



C 



A C* 
C 



\4 


C*~ 




'A 


B 


C 





— r 


C 






(1.49) 
(1.50) 
(1.51) 
(1.52) 



The matrices X that satisfy (|1.45j) - (|1.48j) (namely, the global maximizers and minimizers of the objective 
rank and inertia functions) are not necessarily unique and their expressions were also given in [18 by using 
certain simultaneous decomposition of the three given matrices. Observe that the right-hand sides of (|1.45jl - 
(|1.48|) are represented in analytical forms of the ranks and inertias of the five given block matrices, we can easily 
use them to derive extremal ranks and inertias of some general linear and nonlinear matrix- valued functions. 
In these cases, combining the rank and inertia formulas obtained with the assertions in Lemma 1.1 may yield 
various conclusions on algebraic properties of linear and nonlinear matrix- valued functions. 



2 The extremal ranks and inertias of A — B\XBl subject to B^XB^ = 
A 2 and B 3 XB% = A 3 

We first derive a parametric form for the general common Hermitian solution of the pair of matrix equations in 



Lemma 3.1 (|33j) Let Ai G C^ s , Bi s C'"' x " be given for i = 2, 3, and suppose that each of the two matrix 
equations 



B 2 XB* 2 = A 2 and B 3 XB* 3 = A 3 



(3.1) 



has a solution, i.e., 3$ (At) C 3$(Bi) for i = 2, 3. Then, 
(a) The pair of matrix equations have a common Hermitian solution if and only if 



A 2 B 2 
-A 3 B 3 
B* 2 B* 3 



= 2r 



B 2 
B 3 



(3.2) 



(b) Under (|3.2I) . the general common Hermitian solution of the pair of equations can be written in the following 
parametric form 



X = X + VF B + F B V* + F B2 UF B3 + F B3 U*F B2 , 

where Xo is a special Hermitian common solution to the pair of equations, B 
are arbitrary. 



B 2 
B 3 



(3.3) 

and U, V G C nx ™ 



G 



Substituting (O into A 1 - B X XB\ gives 

Ai - BxXBl = A X - BxXoBl - BiVF B B\ - B 1 F B V*B{ - BxF B2 UF Ba B*x - B\ Fb 3 U* Fb 2 B\ , (3.4) 

which is a matrix-valued function involving two variable matrices V and U. Thus, the constrained matrix-valued 
function in (jl.2l) is equivalently converted to the unconstrained matrix- valued function in (|3.4[) . To find the 
global maximal and minimal ranks and partial inertias of (|3.4[) , we need the following result. 



Lemma 3.2 Let 



4>{X 1 , X 2 ) = A - B 1 X l C 1 - (B 1 X 1 C 1 )* - B 2 X 2 C 2 - {B 2 X 2 C 2 ) 



(3.5) 



where A £ CjJ , Bi € C mxpi and Ci € <C qiXm are given, and Xi S C PiX?i are variable matrices for i = 1, 2, and 
assume that 



3>,{B 2 ) C @(Bx), M{C{) C ^(Bi), ^(C£) C M{B X ). 



(3.6) 



N = 



M 



A B 2 C* C* 2 
Ci 



, N 1 = 



A B 2 C* C* 2 
B* 2 
Ci 



A B 2 C\ C* 
d 
C 2 



A B 1 
Ci 



Mi 



A B 2 C{ 
B* 
Ci 



M 2 



4 C{ c 2 * 

Ci 
C 2 



Then, the global maximal and minimal ranks and partial inertias of 4>{Xx, X 2 ) are given by 

mm{r[A, Bx], r(N), r(Mx), r(M 2 )}, 

2r[A, B 1 ]-2r(M)+2r(iV)+max{s 1) s 2) s 3 , «4 }, 

min{i ± (Mi), i±(M 2 )}, 



max r[0(Xi,X 2 )' 

X i GCPl x 91 , X 2 Gt>2 X 92 

min r[<j)(Xi, X 2 ) 

Xi eCi x 91 , X 2 6C p 2 x «2 



max X 2 ) 

Xi ec^i x 91 , x 2 ec p 2 x 72 

min i± [ , X 2 ) ] = r[ A, B x ] - r(M) + r{N) 

Xi SC p l x 91 , X 2 £C»! x 92 



+ max{z±(M 1 ) - r(iVi), i±(Af 2 ) - r(N 2 ) }, 



(3.7) 
(3.8) 
(3-9) 

(3.10) 



where 



ax = r(Mi) - 2r(2V x ), s 2 = r(M 2 ) - 2r(JV 2 ), 
s 3 = i+(Mi) + i_(M 2 ) - r(JVi) - r(N 2 ), 
s 4 = i+(Mi) + i_(M 2 ) - r(JVi) - r(N 2 ). 

Proof Under (|3.6|) . applying Lemma Tl. 121 to the variable matrix X\ in (|3.5p and simplifying, we obtain 

^4 — B 2 X 2 C 2 — (B 2 X 2 C 2 )* CI 



maxrMXt, X 2 ) ] = min { r[A - B 2 X 2 C 2 - (B 2 X 2 C 2 )*, Bx], 

Xi 



Cx 







min < r[A, Bx], 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 )* CI 
Cx 



mmr[(j)(Xx, X 2 )] = 2r[A - B 2 X 2 C 2 - (B 2 X 2 C 2 )* , Bx 

Xi 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 ) 
Cx 



ci 





maxi±[4>(Xx, X 2 ) ] = i± 

X\ 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 )* CI 
Cx 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 )* Bx 
Cx 



- 2r 
2r[A, Bx]+r 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 )* CI 







2r 



A Bx 
Cx 



(3.11) 



(3.12) 



(3.13) 



xmxxi ± [4>(Xx, X 2 )' 

Xi 



= r[A-B 2 X 2 C 2 -(B 2 X 2 C 2 )*, B ± ]+i± 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 ) 
Cx 



Cx* 




r[A, B 



A — B 2 X 2 C 2 — (B 2 X 2 C 2 )* Bx 
Cx 

A — B 2 X 2 C 2 — (B 2 X 2 C 2 ) 
Cx 



■ i± 



cf 




' A 


Bx' 





— r 








(3.14) 
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Notice that 



a v ( d v \* 

A — -D 2 A 2 2 — \t>2-r^2 ( -'l) L'J 








£>2 


Ci 




Ci 











V>(^ 2 ). 







X 2 [C2, 0] 



C3 



X 2 *[B 2 *,0] 



Applying Lemma 11.111 to this expression gives 



max r[V>(X 2 )l = minlrW), r(M\), r(M 2 ) }, 

X2gC mx "2 

min r[?/;(X 2 )] = 2r(N) + max{ sx, s 2 , s 3 , S4}, 

X 2 EC mx "2 

max i±[^(X 2 )] = min{ z±(Mi), i±(M 2 )}, 

x 2 ec mx P2 

min i±[^(X 2 )] =r(JV) + max{i±(M 1 ) -r(Nx), i±(M 2 ) - r(N 2 ) } , 

X 2gC mXp 2 



(3.15) 

(3.16) 
(3.17) 
(3.18) 
(3.19) 



where 



Sl = r{Mi) - 2r{N x ), s 2 = r(M 2 ) - 2r(N 2 ), 

s 3 = i+(Mi) + i_(M 2 ) - r(2Vi) - r(JV 2 ), s 4 = i_(Mi) + i+(Af 2 ) - r(JV x ) - r(JV 2 ). 

Substituting these results into (|3TTT]) - (pn3|) yields (j3~7| - ([3~TU| . □ 

It is obviously of great importance to be able to give analytical formulas for calculating the global maximal 
and minimal ranks and inertias of the matrix expression in p.5[) under the assumptions in p.7[) . However, it 
is not easy to find the global maximal and minimal ranks and inertias of a general <fr(Xi, X 2 ) as given in (13. 5[) . 
For convenience of representation, we rewrite Q3.4p as 

Ai - BiXBl = A- GxVG 2 - (G 1 VG 2 )* - G 3 UG 4 - (G 3 UG 4 )*, (3.20) 

where 

A = Ax — BiX B*, G\=B\, G 2 = FbB±, G 3 = BiBb 2 , G4 = B^Bf. (3-21) 
It is easy to verify that the above matrices satisfy the conditions 

(G 2 ) C#(Gi), Sl{Gz)Q&(Gx), 8%{G%) C M(Gx), M{G* 2 ) C #(G a ), ^(G 2 ) C ^(Q). (3.22) 
In this case, applying Lemma T3. 21 to (|3.22p yields the main results of this section. 

' 4 and Bi £ C m,x ™ oe given for i = 1, 2, 3, and assume that the pair of matrix 



Theorem 3.3 Let ,4, e ( 

equations 

B2XB2 = A 2 and B3XBI = A 3 
have a common solution X G C§. Also denote the set of all their common Hermitian solutions by 

S = {X G Cg I B 2 XB* = A 2 , Sa-^-BS = A 3 }. 

and Zet 



Pi 



A x B x 







Bi* B 2 B| 



B> 



Ax 





Pi 





-A 2 


B 2 


/; 


B* 2 






Pi 



Ax 


BJ B 3 * 



Ai 







-A 2 




Bi Bi 
B 2 

-A 3 B 3 



B* B* 



PI 







Q.2 = 



Ax 





Bi 


Bi 





-A 2 


B 2 





Bi* 


B 2 * 

















B 3 



Q3 



Bi 

A 3 B 3 


A x 








Bi Bi 



-A, B, 







B| Bj* 





B 2 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



Then, 



(a) The global maximum rank of Ax — BxXB\ subject to (|3.24[) 



max r (Ai - BxXBt) 
xes 

= mini r[Ax, Bx], r{Qx) 



B 2 
B 3 



r(B 2 ) - r(B 3 ), r(P 2 ) - 2r(B 2 ), r(P 3 ) - 2r(B 3 



(3.27) 
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(b) The global minimum rank of A\ — B\XB\ subject to (|3.24|) is 



minr(Ai - B X XB\ ) 
xes 

= 2r[Ai,fl 1 ]-2r(Pi)+2r(Qi)+max{r(P 2 )-2r(Q 2 ), r(P 3 ) - 2r(Q 3 ), u 1; u 2 }, 



(3.28) 



w/iere 

ui = t+ (P 2 ) + i_(P 3 ) - r(Q a ) - r(Q 3 ), u 2 = i_(P 2 ) + z+(P 3 ) - r(Q 2 ) - r(Q 3 ). 
(c) The global maximum partial inertia of A\ — B\XB* subject to Q3.24p is 



maxi±( A x - B X XB\ ) = min {i±(P 2 ) - r(B 2 ), i ± (P 3 ) - r(B 3 )} . 

JC Go 



(3.29) 



(d) TTie global minimum partial inertia of Ai — B\XB\ subject to (|3.24l) is 



mini±(A 1 -B 1 XB$) = r[A 1 , B x ] - r(P x ) + r(Qi) 



+ max{i±(P 2 ) - r(Q 2 ), i±(Ps) - r(Q 3 )}. 



(3.30) 



Proof Under (|3. 221) . we find by Lemma |3~21 that 



maxrfAi -BtXBl) = maxr[A-G 1 VG 2 ~ (GiVG 2 )* - G 3 UG A - {G 3 UGA*} 
xes v,u 



= min < r[ A, G\ ], r 



A G 3 G\ 
G 2 

minr(Ai - B X XB\ ) = minr [A - G\VG 2 - (Gi^G 2 )* - G 3 [/G 4 - (G 3 t/G 4 )* 
xes v,u 





' A 


G 3 " 






l - T 4 






G% 







G4 





} 



2r[A, Gi] -2? 
+ max{s 1 , s 2 , s 3 , s 4 }, 



A Gi 
G 2 



2r 



A G 3 G4 
G 2 



maxi±( A 1 - BiXBl ) = maxi± [A - G X VG 2 - (G X VG 2 )* - G 3 [/G 4 - (G 3 UG 4 Y 
xes v,u 







G 3 " 




' A 








G* 3 





, i± 


G4 





} 



mini±(Ai -BxXBl) = mini± [A- Gi^G 2 - (Gi^G 2 )* - G 3 J7G 4 - (G 3 [/G 4 )*] 
xes vaj 



•[A, Gi 





' A 


Gi" 




' A 


G 3 




r 


G 2 





+ r 


G 2 









■max{ii, t 2 }, 



(3.31) 



(3.32) 
(3.33) 

(3.34) 



where 



si = r 
s 2 = r 

S3 = i+ 
s 4 = i_ 

*i = i± 
h = i± 



' A 


G 3 " 


- 2r 


G* 3 





' A 


Gl 


- 2r 


G 4 






.4 

G* 3 

A 

G* 3 

A 

G* 3 

A 

G* 



G 3 


G 3 


G 3 ' 






j4 G 3 G4 
G^ 

A G 3 G4 
G 4 

i4 G^ 
G 4 

A G% 
G 4 

A G 3 G 4 
G^ 



.4 

G^ 

' A 

G% 



G 3 






G* 
4 



G* 
4 





A G 3 G! 



G 4 
A G 3 

G4 





G* 
4 





A G 3 
G| 



G* 
4 
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Applying (|1.16|l - (|1.18|l and (|1.25[) . and simplifying by [B 2 X Bf, B 3 X B% ] = [A 2 , A 3 ], elementary matrix 
operations and congruence matrix operations, we obtain 



r[A, G x 

A G 3 G4 
G 2 







A G 3 
G* 3 



= r[A 1 -B 1 X B* 1 , B 1 ]=r[A ll B x ], 

A\ — B\X a B\ BiFb 2 BiFb 3 
FbB* 

A\ — BiXqBI 
Bl 



(3.35) 



Ai 
B* 










Si 


B 2 




Si 


£2 








5:* 



Si 



B 3 

BiXaB*' 
B* 







B* 





"(B) - r(B 2 ) - r(B 3 ) 





Bi 


Bi 








BJ 








B* 


B* 





B 2 





-A 2 











B 3 





-A 3 



- r{B) - r(B 2 ) - r(B 3 ) 



-r(B)-r(B 2 )-r(B 3 ) 



= r{Q 1 )-r(B)-r{B 2 )-r{B 3 ), 



' A 


Gi" 




G 2 





= r 



A\ — B\XqB\ 
F B B\ 

r(Pi)-r(S), 



Bi 




Bt 



Bi 






B* 



-r(B) 



(3.36) 



(3.37) 



Ai — B\X a B\ 
Fb 2 B{ 



B\Fb 2 




= i± 



BiX B* 
Bl 




B* 

B1X0B2/2 
i±(P 2 )-r(B 2 ), 



Bi 


B 2 



B1X0B2/2 
B* 




-r(B 2 ) =i± 



Bi 


B 2 

Ai 
B* 






/' 



Bi 


B 2 



- r(B 2 ) 




B 



2 
A 2 



- r(B 2 ) 



(3.38) 



A G 3 G| 
G^ 



Ai — BiXqBI B\Fb 2 B\Fb 
F B2 B{ 

Ai-BiX B* Bi Bi 



Bl 









BX 



B 2 
B 3 



- 2r(B 2 ) - r(B 3 ) 



A\ B\ B\ B\XqB% 
Bl 
B 2 



Bl 

B 3 

Ai Bi Bi 

B* B, 

B 2 

B 3 

r(Q 2 ) - 2r(B 2 ) - r(B 3 



2r(B 2 ) - r(B 3 



r 



2 

A 2 




2r(B 2 ) - r(B 3 



By a similar approach, we can obtain 



A d 
G* A 



i±(P 3 )-r(B 3 ), 



A G 3 G* A 



Ga 







r(Q 3 ) - r(B 2 ) - 2r(B 3 ). 



Substituting ([3~35]) - ([3T40]) into (f3T3T]) - (f3734|) yields ([3^27]) - ([330]) . □ 

Some direct consequences of the previous theorem are given below. 



(3.39) 



(3.40) 



Corollary 3.4 Let Ai € 



and Bi 6 



&e given for i = 1, 2, 3, and suppose that each pair of B\XBl 
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Ai, B 2 XB% — A2 and B3X B3 = A 3 have a common Hermitian solution. Also let S be of the form (|3.23|) . Then, 

r(B 2 )-r(B 3 ), 
2r(B 3 ) 



maxr( A 1 - B x XBl) = min<J r{Bi), r{Qi) - r 



B 2 
B 3 



2r 



B 2 



2r(B 2 ), 2r 



min r{ A x - B X XB\ ) = 2r(Q 1 ) - 2r 



xes 



Bi 
B 2 
B 3 



2r 





'Bi 







-r(B 2 ), 



£1 
B 3 

Bi Bi 
B 2 
B 3 

Bi' 

B 3 



r(B 3 ) , 



xes 





Bi 




Bi 


Bi' 


BxXBZ) = r(Qi)-r 


B 2 


— r 


B 2 







B 3 







B 3 



(3.41) 
(3.42) 
(3.43) 
(3.44) 



where Qi is of the form (I3.26[) . 



Proof Under the given conditions, the ranks and inertias of the block matrices in (|3 . 25[) and (|3.26[) are given 
by 



r(Pi) =r{Bi) +r 



Bi 
B 2 
B 3 



r(P 2 ) = 2r 



Bi 
B 2 



, r(P 3 ) = 2r 



Bi 

B 3 



r{Q 2 ) = r 

Hence (j3~27| -(j3~30 l) reduce to (j3T41"1) -(|3T44 j) . □ 



Bi 


Bi' 




B 2 





+ r 





B 3 





Bi 
B 2 



r(Q 3 



, i±(P2)=r 



Bi Bi 
B 2 
B 3 



Bi 
B 2 



, i±(P 3 )=r 



Bi 
B 3 



Bi 
B 3 



Corollary 3.5 Let Ai € Cg' xm ' and Bi G c^iX" ^ e gi ven for i = 1, 2, 3, and suppose that each pair of the 
triple matrix equations 

BiXB{ = A x , B 2 XB; = A 2 , B 3 XB* = A 3 (3.45) 
have a common Hermitian solution. Then, there exists a Hermitian X such that (|3 .45[) holds if and only if 



Ai 








Bi 


Bi 





-A 2 





B 2 











-A 3 





S3 


Bl 


B* 2 


B* 3 









Bi Bi 
B 2 
B 3 



r[B*i,B* 2 ,B* 3 [ 



(3.46) 



Proof It follows from (pT42|) . □ 

A challenging open problem on the triple matrix equations in (|3.45|) is to give a parametric form for their 
general common Hermitian solution. 

Setting Bi — I n in Theorem 3.3 may yield a group of results on the extremal ranks and inertias of Ai — X 
subject to (|3.24j) . In particular, we have the following consequences. 

Corollary 3.6 Let Ai 6 C^ 1 and Bi £ c m<xn be given for i = 2, 3, and assume that (j3.23|) has a common 
solution. Also let S be of the form (13.24|) . Then, 



(a) The global maximum rank of the solution of (|3.24| is 



maxr(J) = min{ n, s x , s 2 , S3}, 

X tE<S 



(3.47) 



where 



si — 2n - 


- r 


' A 2 



B 2 
A 3 B 3 


— r 


' B 2 


s 2 — 2n- 


V r(A 2 ) 


-2r(B 2 ), 


s 3 = 


= 2n 



r(B 2 ) - r(B 3 ), 
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(b) The global minimum rank of the solution of (|3.24l) is 



min r(X) = 2r 



A 2 B 2 
A 3 B 3 



■max{ti, t 2 , t 3 , i 4 }, 



where 



t,-K4)->r[t | 
t 3 = i+(A 2 ) +i-(A 3 ) -r 

U = i-(A 2 ) + i+(A 3 ) - r 



, h = r(A 3 ) - 2r 



A 2 B 2 

B 3 

A 2 B 2 

B 3 



B 2 
A 3 B 3 

B 2 
A 3 B 3 

B 2 
A 3 B 3 



(c) TTie global maximum partial inertia of the solution of (|3.24[) is 

maxi±pO = min{ri + ! ± (A 2 ) -r{B 2 ), n + i±(A 3 ) ~ r(B 3 )}. 

(d) The global minimum partial inertia of the solution of (|3.24|) is 



min i+(X) = r 



A 2 B 2 
A 3 B 3 



+ max I i±(A 2 ) - r 



A 2 B 2 
B 3 



, *±(40-r 



B 

A 3 B 



In consequence, 
(e) Eq. (|3.23p raas a solution X > if and only if 

A 2 > 0, A 3 ^ 0, <£(A 2 ) = #(fl 2 ), ^(A 3 ) - «(B 3 )- 



(f ) All solutions of (|3.23p satisfy X > if and only if A 2 ^ 0, A 3 ^ and one of 

r(A 2 ) = r(B 2 ) = n, r(A 3 ) = r(S 3 ) = »■ 



(g) i?g. (|3.23p has a solution X < if and only if 

A 2 ^0, A 3 ^0, %{A 2 )=M{B 2 ), M{A 3 )=M{B 3 ). 



(h) All solutions of (|3.23p satisfy X < if and only if A 2 0, A 3 ^ and one o/ 

r(A 2 ) = r{B 2 ) = n, r(A 3 ) = r{B 3 ) = n. 



(i) Eq. (|3.23p has a solution X ^ if and only if 



A 2 > 0, A 3 ^ 0, 



A 2 




c M 



5 2 
A 3 B 3 





-4 :! 



c 



A 2 B 2 
B 



(j) All solutions of (|3.23p satisfy X ^ if and only if A 2 ^ 0, A 3 and one o/ 

r {B 2 ) = n and r(B 3 ) = n. 



(k) Eq. (|3.23p raas a solution X ^ if and only if 



A 2 € 0, A* < 0, St 



A 2 




c ^ 



B 2 

^3 ^3 





M 



c 



3 



A 2 B 2 
B 3 



(1) >U/ solutions of (|3.23p satisfy X ^ if and only if A 2 ^ 0, A 3 ^ and one o/ 

r (B 2 ) = n and r(B 3 ) = n. 
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Proof Set A\ — and B\ = I n in Theorem 13.31 and simplifying, we obtain (a)-(d). Applying Lemma Tl. 5 1 to 
(QLll and (pT39]l . we obtain (e)-(l). □ 

Corollary I3.6f e s )-(1) give a set of analytical characterizations for the existence of definite common solutions 
of the two matrix equations in (|3.23[) by using some rank and range equalities and inequalities. These character- 
izations are simple and easy to understand in comparison with various known conditions (see, e.g., |14 [ I40ll41"] )s 
on the existence of definite common solutions of (|3.23p . 

Rewrite B 2 XB\ = A 2 and B 3 XB% = A 3 as 



B21, B 2 



X\ X 2 

XX X, 



B. 

B% 



21 



= Ao 



. B 31 , B 32 ] 



X\ X 2 

x* 2 x 3 



TDM 

B31 
B 32 



= A 3l 



(3.51) 



where B a G C m ' xni , B i2 G 



2, 3, Xi G 



-H i 



X 2 G C" lX " 2 and X 3 G C^ 2 with m + n 2 = n. We 



next derive the extremal ranks and inertias of the submatrices X\ and X 3 in a Hermitian solution of (|3.5ip . 
Note that X±, X 2 , X 3 in ()3.51j) can be rewritten as 



Xi = PxXPf, x 2 = Pixp;, x 3 = p 2 xp;, 



(3.52) 



where Pi = [I ni) 0] and P 2 — [0, I n2 ]. For convenience, we adopt the following notation for the collections of 
the submatrices X\ and X 3 in (|3.51|) : 



S 1 = {X, = P,XP* 1 B 2 XB* = A 2 , b 3 xb; 

s 3 = {x 3 = p 2 xp; 1 b 2 xb; = a 2 , b 3 xb. 



3 



A3, x 
A 3 , X 



x*}, 

X*}. 



(3.53) 
(3.54) 



The global maximal and minimal ranks and partial inertias of the submatrices X\ and X 3 in (|3.5ip can easily 
be derived from Theorem 13.31 The details are omitted. 

If each of the triple matrix equations in (|1.8j) is not consistent, people may alternatively seek its common 
approximation solutions under various given optimal criteria. One of the most useful approximation solutions 
of BXB* = A is the well-known least-squares Hermitian solution, which is defined to be a Hermitian matrix X 
that minimizes the objective function: 



|| A- BXB* || 2 = ti[( A - BXB* )( A- BXB*)*}. 
The normal equation corresponding to the norm minimization problem is given by 

B*BXB*B = B*AB. 



(3.55) 



(3.56) 



This equation is always consistent. Concerning the common least-squares Hermitian solution of (jl.8p . we have 
the following result. 

Corollary 3.7 Let Ai G CJJ* and Bi G C miXn be given for i = 1, 2, 3. Then, triple matrix equations have a 
common least-squares Hermitian solution, namely, there exists an X G C^ x ™ such that 



Ai - BiXB* 



I, 2, 3, 



(3.57) 



if and only if 



B*A t Bi 


B*B. t 

B\A X B X 



B\B X 





B* Bj 


-B*A 2 B 2 


BXBo 



BtBi 
B* Bj 




2r 



Bi 
B, 






-B* 3 A 3 B 3 
B%B 3 



B\B X B\B X 
B*B 2 
B*B, 



> 2 ^ 2 





3-^3 





i, j = 


1, 2, 


3, 








Bi 


Bx~ 




~B{ 


= r 


B 2 





+ r 


B 2 







B 3 




B 3 



(3.58) 



(3.59) 



Proof It follows from Lemma 13. li Corollary 13.51 and (|3.56l) . 



4 The extremal ranks and inertias of A\ — B\XB\ subject to the 
Hermitian solutions of B^X = 

Also B4X = A4 in (|1.3p is not given in symmetric pattern, it may have a Hermitian solution, as shown in 
Theorem 1 1.101 So that the global extremal ranks and inertias of A\ — B\XB\ subject to the Hermitian solution 
or nonnegative definite solution of B^X = A4 can also be derived. 
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Theorem 4.1 Assume that the matrix equation B 4 X = A 4 in (|1.3[) has a Hermitian solution, i.e., ^%(A 4 ) C 
M(Bi) and A 4 B\ = B 4 A* 4 , and let 



S = {X e C£ I B 4 X = A 4 }, M 



A^B^ B^ 



N 



Ax Bi 
B\ 



B 4 -A4BI 



Then, 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 

to 

In consequences, 

(a) = A4 /ias a solution X G smc/i i/iai Ai —B\XB\ is nonsingular if and only ifr(M) — r(B 4 ) + mi. 

(b) Ai — B^XB\ is nonsingular for all Hermitian solution of B 4 X — A 4 if and only if 2r{M) = r(N) + mi- 

(c) The pair of matrix equations B\XB\ = A\ and B 4 X — A 4 have a common Hermitian solution if and 



maxr( A\ 
xes 


- BiXB* ) 


= r(M) - 


r(B 4 ), 


min r( A\ 

xes 


- BiXB* ) 


= 2r(M) 


-r(N), 


maxi+( Ai 
xes 


- BiXB* ) 


= i±(JV)- 


- r(B 4 ) 


min i-i-( A\ 

xes 


- B Y XB\ ) 


= r(M) - 


t T (N). 



only if M 



A 1 
A 4 B\ 



c m 



(d) B\XB1 = Ai holds for all Hermitian solutions of B 4 X = A 4 if and only if r(M) = r{B 4 ). 

(c) B 4 X = A 4 has a solution X G Cg such that Ai - BiXB* > {A x - B X XB* X < 0) if and only if 
i+{N) = r{B 4 ) + rrn ( i_[N) = r(B 4 ) + m x ). 

(f) Ai — B X XB* > ( Ai — BiXB\ < 0) holds for all Hermitian solutions of B 4 X — A4 if and only if 
r(M) = i-{N) + mi (r(M) = i + (N) + m 4 ). 

(g) B 4 X = A 4 has a solution X G such that Ai - B X XB\ > ( A x - B X XB\ ^ ) if and only if 
r(M) = i+{N) (r(M) = i_(iV)) . 

(f) Ai — B 4 XB\ ^ ( Ai — BiXB* ^ 0) holds for all Hermitian solutions of B 4 X — A 4 if and only if 
i_{N) = r{B 4 ) {i+(N)=r{B 4 )). 

Proof. From Lemma 11. 1 Of a), the general Hermitian solution of B 4 X = A 4 can be written as 

X = B\A 4 + (B\A 4 )* - B\A 4 B\B 4 + F Bi WF Bl , (4.6) 

where W G is arbitrary. Substituting (|4.5[) into A\ - BiXB* gives 

Ai - B X XB\ =P- BiF Bi WF Bi Bl, (4.7) 

where G = A x - B 4 b\A 4 B\ - B 1 {BlA 4 )*Bl + BiB\A 4 B\B 4 B\. Applying (fi~53V(fT35]) to gU) yields 

r[G, BiF Bi ], 



maxr( A 4 
xes 


- BiXB\ ) 


= max r(G 

wee™ 


-BiF Bi WF Bi Bl 


min r( Ai 

xes 


- BiXB* ) 


= min r( G 

weCft 


-BxF Bi WF Bi Bl 


maxi±( Ai 
xes 


- BiXB* ) 


= max r(G 

wec^ 


-BiF Bi WF Bi B* 


min i±( Ai 

xes 


- BiXB* ) 


= min r( G 


-BiF Bi WF Bi B* 



G BiF Bi 
F Ba B$ 



G B 4 F Bi 
F Bl B* 



G BiF Bi 
F Bi B* 



(4.8) 
(4.9) 
(4.10) 
(4.11) 



It is easy to verify that under B 4 B\A 4 = A 4 , the equality B 4 {b\A 4 )* = B 4 A%(B\)* = A 4 Bl(B\)* = A 4 B\B 4 
holds. In this case, applying It is easy to verify by (1-17) and (1-25) to (|4.8p - (|4.11[) and simplifying by elementary 
matrix operations and congruence matrix operations, we obtain 

~A X - BiB\A 4 B{ - Bi(B\A 4 )* B{ + B 4 b\A 4 b\b 4 B\ B x 

B 4 

Ai Bi 
A 4 B\ + B 4 {B\A 4 )*Bl - A 4 B\B 4 B* Y B 4 

Ai Bi 



r[G, B\F Bi ] = r 



r(B 4 



r{B 4 ) 



A^B* B4 



- r(B 4 ) = r{M) - r(B 4 ), 



(4.12) 
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G BiFb 4 
F Bi Bl 

Ax - BxB\A a BI - Bi(b\A a )* B{ + BxB\A 4 B\B 4 B\ B x 



= i± 



Bl 


A x Si 
B{ 
^A 4 B^ + ^B 4 Bx{B 4 A 4 )* Bi — ^A 4 b\b 4 B\ B4 

r(B 4 ) 







Bl 
B A 



-r(B 4 



iB\B\A 4 B\ + ~^B\A\ — }^B\B\A 4 B\ 
Bl 



B{ 
^A 4 B^ B 4 



Bx ^B\A\ 



Bl 



- r(fl 4 ) = i± 



BJ 



Bi -A4BI 



r(B 4 ) = i±{N)-r{B 4 ). 



(4.13) 



Substituting (|4TT2|) and (|4~T3)) into (j4T8| - (|4lTt yields (|4T2 )) -(|4T5 l) , Applying Lemma Q to ([4T2 ]) -([43 ]l yields 
(a)-(f). □ 



Theorem 4.2 Assume that the matrix equation B 4 X — A 4 in <| 1 .3|) has a nonnegative definite solution, i.e., 
M(A 4 ) C ^(S 4 ), A 4 .B| ^ and r(A 4j B|) = r{A 4 ), and let 



5 = {0CXe 



Then, 



A 4 X = B 4 }, Mi 



Ai Si 



Af 2 



Ai B X A| 
A 4 B^ A 4 B* 4 



N 



Ax Bx 
Bl 
B 4 





B| 
-A 4 B 



4-D4 
(4.14) 



maxr( Ai 
xes 


- BxXB* ) 


= r(Mx) - 


r(B 4 ), 


min r( Ai 


- BxXB{ ) 


= r(Mi) + 


i_(M a ) 


max i + ( Ai 


- BxXB'l ) 


= i+(M 2 ) 




min i + ( Ax 

xes 


- BxXB\ ) 


= r(Mi) - 


i_(JV), 


maxi_( Ax 
xes v 


- BxXB{ ) 


= i_(JV)- 


r(At), 


min z_( Ax 

xes 


- BxXB'l ) 


= i_(M a ). 





(4.15) 

_(iV), (4.16) 

(4.17) 
(4.18) 
(4.19) 
(4.20) 

A to 

In consequences, 

(a) B 4 X = A 4 has a nonnegative definite solution such that Ax — BxXB\ is nonsingular if and only if 
r(Mx) = r{B 4 )+mx. 

(b) Ax — BxXB\ is nonsingular for all nonnegative definite solution of B 4 X = A 4 if and only if r(Mx) + 
i-(M 2 ) = i-(N) + m x . 

(c) The pair of matrix equations BxXB* = Ax and B 4 X = A 4 have a common nonnegative definite solution 
if and only if 'r(Mi) + i-(M 2 ) = i-(N). 

(d) BxXB\ — Ax holds for all nonnegative definite solutions of B 4 X = A 4 if and only if r(M) — r(B 4 ). 

(e) B 4 X = A 4 has a solution X € C\\ such that Ax — BxXB[* > if and only if i + (M 2 ) = r{A 4 ) + m±. 

(f) Ax — BxXBi > holds for all Hermitian solutions of B 4 X = A 4 if and only if r(Mx) = i-(N) + mi. 

(g) B 4 X = A 4 has a solution X £ Cg such that Ax — BxXB\ < if and only if i-(N) = r(A 4 ) + mi. 

(h) Ax — BxXB\ < holds for all Hermitian solutions of B 4 X = A 4 if and only if i-(M 2 ) = Tn x - 

(i) B 4 X = A 4 has a solution X € C\\ such that Ax — BxXB\ ^ if and only if M 2 ^ 0. 

(j) Ax — BxXB\ ^ holds for all Hermitian solutions of B 4 X = A 4 if and only if i^(N) = r(A 4 ). 
(k) B 4 X = A 4 has a solution X G Cg such that Ax - BxXB{ sC if and only if >(Mi) = i-(N). 
(1) Ax — BxXB\ ^ holds for all Hermitian solutions of B 4 X = A 4 if and only if i^{M 2 ) — r(A 4 ). 
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Proof. From Lemma ll.lOf b). the general nonnegative definite solution of B 4 X = A 4 can be written as 

X = Al{A 4 B*^A 4 + F Bi WF Bl , (4.21) 
where W e Cg is arbitrary. Substituting gl| into A 1 - B X XB\ gives 

A x - B X XB\ = G — BxFmWFinBl, (4.22) 
where G = A x - BxA\(A 4 B\)1 A 4 B\. Applying (|i~4TJ)) - ([P2"|) to (|4~22l yields 



maxrfii - B X XB*) 

XGS 



min r( Ax - BxXB* ) 



maxi + (Ax - BxXB* ) 



mini + (Ai - BxXB* ) 



maxi_(Ai -BxXB*) 
xes 

mmi±(Ax — BxXBf ) 

X<ES 



max r(G- B\ Fb 4 WF Bl B* )=r[G, BxF Bi }, 

min r(G-BxF B4 WF Bi Bt)=i-(G)+r[G, B x F Bi ]-i. 
o^wecv, 



(4.23) 



G 
F Bi Bt 



BxF Bi 



max r{G- Bx F Bl WF Bi £?* ) = i+{G), 
o^wecv, 



max i+( G - B\ F Bi WF Bi B\ )=r[G, B X F B , ] - z_ 



G BxF Bi 
F Bi B* 



max i_( G - BxF Bl WF Bi Bl ) = L 



G B\F Bi 
F Bi B\ 



min i_( G - BxF Bi WF Bi B* ) = i_(G). 



It is easy to verify by (1.17), (1.23) and (1.25) that 
r[G, BxF Bi ]=r 



Ax - BxAKMBtV MBl Bx 
B 4 



- r(B 4 ) = 



Ax Bx 



-r(B 4 ), 



(4.24) 
(4.25) 

(4.26) 

(4.27) 
(4.28) 



(4.29) 



G BxF Bi 
F Bi B*x 



= i± 



Ax - BxA%(AiBty A 4 B{ Bx 



B{ 


Ax Bx \BxA\ 



B\ 
^A 4 B* B 4 



Bl 



B\ 
B 4 

- r(B 4 ) = i± 



- r-(B 4 ) 



Ai Bx 
Bx* 







i ± (G) = »±[Ai - BtAJ^SiJt^B* ] = i ± 



P 4 



«±(^ 4 B 4 *) 



(4.30) 



(4.31) 



Substituting (j4T2^|) - (jOT|) into (|4T2l|) - (|4T2"gj) yields P35l - (|4T2U1) . Applying Lemma O to g^-g^]) yields 
(a)-(l). □ 

Corollary 4.3 Assume that the matrix equation in Lemma 1.10 has a Hermitian solution, P € Cg, and let 
S = { X G Cg I AX = B }. Then, 



maxr(X-P) = r(B-AP) -r(A)+n, 

X £0 

min r( X - P ) = 2r( B - AP ) - r{ BA* - APA* ), 

X OzS 

maxi±(X - P ) = i±(BA* - APA*) -r(A)+n, 
xes 

mmi±(X - P) = r(B - AP ) — i T (BA* - APA* ). 

X £tS 



(4.32) 
(4.33) 
(4.34) 
(4.35) 



In consequence, 

(a) There exists an X e S such that X — P is nonsingular if and only if 

M(AP-B) =M{A). 

(b) X — P is nonsingular for all X £ S if and only if 

2r( B - AP ) = r( BA* - APA* ) + n. 

(c) There exists an X 6 S such that X > P (X < P) holds if and only if 

@{ BA* - APA* ) = M(A) and BA* ^ APA* ( ^( BA* - APA* ) = M(A) and BA* ^ APA* 
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(d) X > P (X < P) holds for allX eS if and only if 

r(B-AP) = n and BA* > APA* (r(B - AP) =n and AB* < APA* ) . 

(c) There exists an X G S such that X ^ P (X ^ P) holds if and only if 

&(B—AP) = &(BA*—APA* ) and BA* > APA* (&(B — AP ) = BA* — APA* ) and BA* < APA* ). 

(f) X ^ P (X P) holds for allX eS if and only if 

BA* ^ APA* and r(A) = n (BA* < APA* and r(A) = n) . 

Corollary 4.4 Assume that the matrix equation in Lemma 1.10 has a Hermitian solution X ^ 0, and let 
O^Peq. Also, define 



S = {OsC X gC£| AX = B}, M = 



BA* B 
B* P 



(4.36) 



Then, 



maxrf X - 
xes 

min r( X — 

xes 


P) 
P) 


= r(B-AP) -r(A)+n, 

= i-(M) + r(B-AP)-i+{BA* -APA*), 


(4.37) 
(4.38) 


maxi+( X — 
xes 


P) 


= i+(BA* - APA* ) - r(A) + n, 


(4.39) 


min i+( X — 

xes 


P) 


= i-(M), 


(4.40) 


maxi_( X — 
xes 


P) 


= i+(M)-r(B), 


(4.41) 


min i_( X — 

xes 


P) 


= r(B - AP) - i + (BA* - APA* ). 


(4.42) 



In consequence, 



(a 
(b 
(c 
(d 
(e 
(f 
(g 

(n 

(i 
(j 



There exists an X G S such that X — P is nonsingular if and only if B — AP ) = $%(A). 
X - P is nonsingular for all X e S if and only if i_(M) +r(B - AP) = i + (BA* - APA* ) + n. 
There exists an X G S such that X > P holds if and only if3!( B A* - APA* ) = M(A) and BA* ^ APA*. 
X > P holds for all X G S if and only ifi-(M) = r(A). 

There exists an X G S such that X < P holds if and only ifi-(M) — r(B) + n. 
X < P holds for allX eS if and only if r(B - AP) = n and BA* ^ APA*. 

There exists an X G S such that X ^ P if and only if M(B-AP) = 0Z( BA* - APA* ) andBA* ^ APA*. 

X^P holds for allX ES if and only ifi-(M) = r(B). 

There exists an X G S such that X < P if and only if M ^ 0. 

X < P holds for allX eS if and only if i + (BA* - APA* )=n- r(A). 



Corollary 4.5 Assume that the matrix equation in Lemma 1.10 has a Hermitian solution. Then, 

Ax^xec- r{X) =n + r{B) ~ r(A) ' (4 ' 43) 

4 v ^-cr» r{X) = 2r(B) ~ r(AB * } ' (4 ' 44) 

max i ± (X) =n + i ± {AB*) - r(A), (4.45) 

/iA— .fcf, A fc'L-j^ 

A r r.cr. 1±{X) = r{B) ~ 1 t(AB*). (4.46) 

Hence, 

(a) AX = B has a nonsingular Hermitian solution if and only ifr(A) = r(B). 
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(b) AX = B has a solution X > (X < 0) if and only if AB* ^ and r(AB*) = r(A) (AB* ^ and 
r{AB*) = r(A)). 

(c) AX = B has a solution X > (X ^ 0) if and only if AB* > and r(AB*) = r(P) (AB* sC and 
r(AB*) = r{B)). 

(d) The rank of the Hermitian solution of AX = B is invariant the positive index of inertia of the Hermitian 
solution of AX — B is invariant the negative index of inertia of the Hermitian solution of AX = B is 
invariant <^ r(AB*) — r(A) + r(B) — n. 

Finally, we rewrite the matrix equation AX = B as 



'A x , A 2 



Xi X2 
X* x 3 



= [Bi, P 2 ] , 



(4.47) 



where A, G C mx '\ Bi G C mx ™*, X x G Cg 1 , X 2 G C" lX ™ 2 , X 3 G Cg 2 for i = 1, 2 and m +n 2 = n. Note that 
the unknown submatrices in (|4.47l) can be written as 



X\ — P\X P^ , X 2 — PiXP 2 , X3 — P 2 XP^ 



2 > 



(4.48) 



where Pi = [I ni , 0] and P2 = [0, I n2 ]• We next find the extremal ranks and inertias of the submatrices X\ 
and X3 in a Hermitian solution of (|4.47j) . For convenience, let 

Si = {Xi G C^ 1 I Xi = PiXP*, AA = P, X G Cg}, (4.49) 
S 3 = {A 3 G C£ 2 I X 3 = P 2 XP*, AX = B, A 6 C£}. (4.50) 

Applying Theorem 14. II to (|4.49[) and (|4.50p gives the following results. The details of the proof are omitted. 

Theorem 4.6 Assume that matrix equation in (|4.17p has a Hermitian solution, and let S\ and S3 be of the 
forms in (|4.49[) and (|4.50[) . Then, the global maximal and minimal ranks and inertias of the Hermitian matrices 
in Si and S3 are given by 



max r(Xi) = m +r[A 2 , Pi] -r(A), 

Xi G<Si 



min r(X-\ 



= 2r[A 2 , Bi]-r 



AB* A 2 



An 



max i± ( X\ ) = m + i± 



AB* Ao 



A% 











-r(A), 



min i±(Xi) =r[A 2 , Bi]-i T 



AB* A 9 



A* 







(4.51) 
(4.52) 
(4.53) 
(4.54) 



and 



max r(X 3 ) = n 2 +r[Ai, B 2 ] - r(A), 
X3ES3 



min r(X 3 )=2r\Ai,B 2 ]-r 
X 3 eS 3 



max i± ( X 3 ) = 112 + i± 
X3ES3 



min i±(X 3 ) =r[Ai, B 2 ]-i T 

X3 G03 



AP* Ai 

A* 

AP* Ai 
A^ 

AB* Ai 

A* 



-r(A), 



(4.55) 
(4.56) 
(4.57) 
(4.58) 



Applying Lemma Tl. 5 1 to (|4.51l) - (|4.54p . we easily obtain the following algebraic properties of the submatrix 
Xi in (14A71) . 

Corollary 4.7 Assume that matrix equation in (|4.4T[) /ias a Hermitian solution. Then, 
(a) (|4.47p has a Hermitian solution in which X\ is nonsingular if and only if r[A 2 , B\ ] = r(A). 



(b) Xi is nonsingular in all Hermitian solutions of (|4.47|) if and only if r 



AB* A 2 
A 



2 



= 2r[A 2 , Bi] 



ni. 



(c) (|4.47p has a Hermitian solution in which X\ > (X\ < 0) if and only if 



AB* A 2 
A* 2 



= r(A) 



AB* An 



A* 







= r(A) 
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(d) X\ > (X\ < 0) in all Hermitian solutions of (|4.47[) if and only if 



AB* A 9 



A% 







r[A 2 , Bi]-m 







AB* 


A 2 




(<♦ 


A*. 






r[A 2 , Bi]-m 



(e) (|4.47p has a Hermitian solution in which X\ (X\ ^ 0) if and only if 



AB* A 2 



A* 







r[A 2 , S x ] 





AB* 


A 2 


(«- 


A* 






r[A 2 , Si] 



(f) Xi (Xi ^ 0) m all Hermitian solutions of (|4.47[) if and oni?/ if 

AS" 



4* 



A 2 




r(A) — ni 







AS* 


A 2 


ni 


(* 


. A* 2 






= r(A) - n-j 



(g) (|4.47p has a Hermitian solution in which X\ — if and only if ffl(B\) C 3£(A 2 ). 

(h) X\ = in all Hermitian solutions of (|4.47|) if and onZy if r[ A 2 , B\ ] = r(A) — n\. 

(i) T/ie ranfc of X\ in the Hermitian solution of (|4.47p is invariant the positive index of inertia of X\ 
in the Hermitian solution of (|4.47p is invariant O the negative index of inertia of X\ in the Hermitian 

AB* Ao r . ^i / 

- =r{A 2 ,B 1 )+r{A)-n x . 



solution of (|4.47p is invariant r 







5 Conclusions 

In this paper, we studied the problems of maximizing and minimizing the rank and inertia of the constrained 
matrix expression in (11.21) and (jl.3[) . and obtained many symbolic formulas for calculating the extremal ranks 
and inertias of (11.21) and (|1.3|) by using pure algebraic operations of matrices and their generalized inverses. 
As direct applications, we gave necessary and sufficient conditions for the existence of X satisfying the triple 
matrix equations in (|1.2p and (|1.3p . as well as some matrix inequalities. Although the problems of maximizing 
and minimizing ranks and inertias of matrices are generally regarded as NP-hard, the results presented in this 
previous sections as well as the papers [131 [TSl HH [T7J [TSJ 1301 (S3 1321 IM1 13H] show that many closed- form formulas 
for calculating global extremal ranks inertias of some simpler matrix expressions can be established symbolically 
by using some pure algebraic operations of matrices, while these explicit formulas can be used to solve many 
fundamental problems in matrix theory, as mentioned in the beginning of this paper. All the results obtained 
in these papers are brand-new, but easy to understand within the scope of elementary linear algebra. This 
series of fruitful researches show that for many basic or classic problems like solvability of matrix equations and 
matrix inequalities, we are still able to establish a variety of innovative results by some new methods. 

Motivated by the fruitful results and the analytical methods used in this paper, we mention some research 
problems for further consideration: 

(a) A challenging task is to give the closed-form for the general common solution of B 2 XB 2 = A 2 and 
S3XS3 = A3 that satisfies X > (< 0, ^ 0, ^ 0), which is equivalent to solving the inequalities 

X Q + VF B + F B V* + F B2 UF B3 + F B3 U*F B2 > (< 0, > 0, < 0). (5.1) 

Moreover, give the extremal rank and partial inertia of A\ — B\XB\ subject to B 2 XB 2 — A 2 and 
B 3 XB* = A 3 and ^ 0. 

(b) Give the extremal ranks and inertias of the LHMF A\—B\XB\ subject to the common Hermitian solution 
of the k — 1 consistent linear matrix equations 

[B 2 XB*,...,B k XB* k ] = [A 2 ,...,A fc ], 
and to establish necessary and sufficient condition for the set of matrix equations 

[ BiXBl , . . . , BkXBl ] = [ Ai , . . . , Afe ] 
to have a common Hermitian solution, as well as a common nonnegative definite solution. 

(c) Give the extremal rank and partial inertia of Ai — B\XB\ subject to a linear matrix inequality B 2 XB 2 ^ 
A 2 . In such a case, it is necessary to first give analytical expression for the general Hermitian solution of 
B 2 XB* 2 ^ A 2 . 
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(d) Give the extremal ranks and inertias of A\ — B\XB\ subject to B 2 X = A 2 and X ^ 0. 

Since linear algebra is a successful theory with essential applications in most scientific fields, the methods 
and results in matrix theory are prototypes of many concepts and content in other advanced branches of 
mathematics. In particular, matrix equations and matrix inequalities in the Lowner partial ordering, as well as 
generalized inverses of matrices were sufficiently extended to their counterparts for operators in a Hilbert space, 
or elements in a ring with involution, and their algebraic properties were extensively studied in the literature. 
In most cases, the conclusions on the complex matrices and their counterparts in general algebraic settings are 
analogous. Also, note that the results in this paper are derived from ordinary algebraic operations of the given 
matrices and their generalized inverses. Hence, it is no doubt that most of the conclusions in this paper can 
trivially be extended to the corresponding equations and inequalities for linear operators on a Hilbert space or 
elements in a ring with involution. 
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